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New modes of the magnetic flux penetration in SQUIDs with ferromagnetic junctions are pre- 
dicted theoretically. Characteristic multinode anharmonicity of the current-phase relation implies 
the coexistence of stable and metastable and tt states for a large interval of strength of the fer- 
romagnetic barrier influence at low temperature. As a consequence, the coexistence of integer and 
half-integer fluxoid conflgurations appears in both rf and dc SQUIDs, and generates two flux jumps 
per one external flux quantum. In dc SQUIDs, new flux jumps are manifested as two dips in the 
critical current dependence on the external magnetic flux. 



I. INTRODUCTION 



A number ofj-iPjOSsibilities to obtain negative coupling between two superconductors in a Josephson junction have 
been proposed.Era The negative sign of the Josephson current corresponds to the coupling energy minimum at the 
phase difference tt between the superconducting electrodes. In superconducting rings containing an odd number 
of such TT shifts, the flow of spontaneous supercurrent may lead to anomalous flux quantization. This has been 
confirmed in experiments with high- T_, superconductors, where the negative coupling is due to d-wave symmetry of 
the superconducting order parameter.cl Also, reversing the direction of the supercurrent appears to be observed in 
out-of-equilibrium mesoscopic superconducting structurcs.u Recently, a metastable tt state has been observed in '^He 
weak links 

The possibilityp-of tt coupling through a barrier containing paramagnetic impurities has been proposed first by 
Brdaevskii et al.,EJ and in superconductor/ferromagnet /superconductor (S/F/S) metallic weak links by Jiuzdju at 
al.El Subseaueatly, the progress in fabrication of S/F multilayers has motivated considerable experimentalSll^ -and. 
theoreticallljLJ investigations of tt coupling. Theory of proximity coupled SylE-|Superlattices in the dirty limit,E£rE£l 
and of atomic-scale S/F superlattices with tunneling between S and F layers ,E2jl3 is well established. A clear evidence 
of TT coupling, which has been sought in several S/F systems from the. oscillations of the critical temperature, predicted 
in Refs. |l^ and now appears to be established unambiguously.EjCJ Recently, another evidence of tt state is found 
in the nonmonotonic variation of the critical current in S/F/S junctions. Lja r— ■ 

The modified Andreev reflection at S/F iatesfaces has been explored theoretically, both for s-waveE3 and d- 
wave superconductors|23 and experimentally.E3'c3 The Josephson current in s-waye, piperconductor/ferromagnetic 
insulator/superconductor (S/FI/S) junctions . h a s .been calculated in several papers,E3"E3 whereas more recent works 
deal with d-wave S/F/S and F/yS j unctions and with the general case of unconventional superconductor junc- 
tions with ferromagnetic barrier .lj 

At low temperature the Andreev reflection is dominant, and. the ferromagnetic barrier may induce a characteristic 
anharmonic current-phasjs relation with new nodes in /(0)E3OE3 The similar type of anharmonicity may appear in d- 
wave pinhole junctions,E3 or_iri S/I/S junctionsJiO as found experimentally in the grain-boundary Josephson junctions 
in high-Tc superconductors.EEl 

In this paper we present new consequences of the modified Andreev reflection in S/F/S junctions. We show that 
coexisting stable and metastable and tt states appear in the vicinity of a critical strength of the ferromagnetic barrier 
influence, Sec. II. This coexistence leads to new modes of the magnetic flux penetration in superconducting loops 
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interrupted by a single ferromagnetic junction (rf SQUID), or two junctions (dc SQUID), Sec. Ill and IV, respectively. 
Section V contains a brief summary and discussion. 



II. CURRENT-PHASE RELATION 



We consider an S/F/S junction with a ferromagnetic metal barrier of thickness d and constant exchange energy 
h. Both superconducting and ferromagnetic metals are assumed clean, with same dispersion relations and with same 
Fermi energies. The pair potential variation is taken as a step function: A = in the barrier, and A — Aq exp(±i0/2) 
in superconducting electrodes, where (f) is the phase difference at the junction. _ 

The Josephson current density is calculated within the quasiclassical theory of superconductivity^] 

j(r) = -2ie7rN{0)kBTj2 (^o^4^) ' ^^'^^ 

where vq is the Fermi velocity, (• • ■) denotes the angular averaging over the Fermi surface, and N{0) is the density 
of states at the Fermi level. The quasiclassical normal Green's functions ga-{r,VoM,i), where a =ti ii and hun — 
TTkBT{2n + 1), are calculated from the Eilenberger equations in the previous paper .e3 

For s-wave pairing in superconducting electrodes and transparent interfaces between two metals, the normal Green's 
function in the barrier is 

_ w„cosx/2 + iV^2 _^ Ag sinx/2 2) 



where 



for a =t, i respectively, and 



y^ul + Ag COS x/2 + ic^n sin x/2 ' 



Z 2iijj„d 

X^^T T — (2.3) 

cos (p Twq cos tf 



2hd , , 

Z = — (2.4) 



is the parameter measuring the ferromagnetic barrier influence, being the angle between the direction of vq and the 
axis perjpendicular to the barrier. For Z = 0, the Green's function, Eq. (2.2), has the well known form for S/N/S 
junctionEa 

For spherical Fermi surface, the Josephson current through the barrier of area 5, I ~ jS, is 

where u — 1/ cos 93, and the normal resistance is given by R^^ = e'^voN{0)S. The temperature dependence of 
Aq = Aq (T) is approximated by 

Ao(T) = Ao(0) tanh (l.74^T,/T - 1^ , (2.6) 

where Tc is the critical temperature of superconducting electrodes. 

The current-phase relation /(</>) is calculated numerically as a function of (j), T and of parameters measuring the 
influence of the ferromagnetic barrier, Z and d. The results, illustrated in Figs. 1 and 2 for fixed thickness of the 
barrier, d = 0.25 ^o(O), where ^o(O) = Tivq/tt Aq (0) is the BCS coherence length at zero temperature, remain practically 
unchanged for any barrier thickness d ~ 0.1 ^o(O) at fixed Z. 

The ground state phase difference (f>gs is determined from the minimum of the S/F/S junction energy 

^0) = !^ r^(0')d0', (2.7) 

where $0 is the flux quantum. Fig. 3. 

At high temperature, T/Tc > 0.9, the current-phase relation is quite close to 



2 



/ = ±/c sin 0, 



(2.8) 



where ± signs correspond to and vr junction, respectively. The amphtude of osciUations in /(0) decreases when Z 
increases from zero to the critical value Zc ~ 1. For Z > Z^^ the supercurrcnt flow in the opposite direction (Fig. 
1(a)), and /c starts to increase. This is the signature of the tt junction with ground state phase difference (/)gs — tt, 
Fig. 3(a). 

At low temperature, for clean and transparent nor ma] m etal barrier, the Andreev reflection is dominant, and 
deviations from the sinusoidal /((/)) law are pronounced.L^cll In the S/F/S case, due to the influence of the exchange 
interaction in the barrier, the current-phase relation is even more anharmonic. Positive for small Z, where state is 
stable, the sign of /(<^) for — > is negative for large Z, where tt state is stable, and multinode oscillations in /(0) 
appear in the vicinity of Z^, Fig. 1(b). In the latter case, there are stationary states /(0) = 0, /'((/i) > both for 
(/> = (27r) and <j> = ir. One of these states is stable, and the other is metastable. For example, at T/Tc = 0.1 for 
Z = 0.9 the metastable state is = tt, = being the stable state, and vice versa for Z — 1.1, Fig. 3(b). Zero 
and TT states are equally stable at Z = Zc- While Zc is practically temperature and barrier thickness independent, 
the interval of the coexistence of and n states, Z^- < Z < Z^-i-, rapidly decreases with increasing temperature. For 
d/£,o{0) — 0.25, Zc = 1.03 and the coexistence interval is 0.1 < Z < 1.5 at T = 0, while Zc = 1.05 and the coexistence 
interval is 1.02 < Z < 1.08 at T/T^ = 0.9. 

The appearence of the n junction for 1 < Z < 4 is also manifested in the nonmonotonous dependence of the critical 
current Ic — max|/((/<)| on Z at given T, Fig. 2. For Z > A, (f)gs = again, andj^car Z ^ 4 metastable states would 
appear as well. However, large Z would induce the decoupling of S electrodes. lj For fixed h, such a nonmonotonic 
variation [©f the critical current as a function of the barrier thickness d has beeij— observed in Nb/Al/Gd/Al/Nb 
junctionsBil The same effect is also predicted for the atomic-scale S/F superlattices£3 

Recently measured I{4>) in ^He weak linksQ are of the same form as the curves in Fig. 1(b) for Z < 1, corresponding 
to the metastable tt state. In d-wave S/F/S junctions, stable and metastable states may also coexist due to the 
exchange interaction. However, in this case the ground state phase differences can varii-Jxom to tt, depending not 
only on Z and T, but also on the crystal orientation of the superconducting electrodesHja 



III. SINGLE-JUNCTION LOOP 



When S/F/S junction is a part of a superconducting loop (rf SQUID), the ground state is determined by minimizing 
the SQUID free energy 



W^(0) = 



$0 

2tt 



/(0')d</>' + ^(0-</.e)' 

47rL 



(3.1) 



where (pc — 27r$c/'I'o represents the normalized external magnetic flux through the loop, and L is the inductance.0 
Therefore, the total magnetic flux through the loop $ = (<I'o/27r)0 is related to the external flux $c by 



L/($) + $ - $o = 0. 



(3.2) 



For <f>c — 0, and for Z sufficiently large (tt junction), there is a flow of spontaneous supercurrent,Lj with nonmonotonic 
dependence on the normalized inductance 



I 



2tt 



-LL 



(3.3) 



Note that when n < <i>c/$o <n+l, where n = ±1, ±2, $/$o ^ ^/^o + n. 

From Eqs. ( ^.5| ) and (3.2) $($c) is calculated numerically. The results are illustrated in Figs. 4 and 5 for several 
values of I and Z, for d = 0.25 ^o(O). 

At high temperature, T/Tc > 0.9, for ^ > 1 and Z < Zc the double-well SQUID energy gives stable and metastable 
solutions of Eq. (3.2), equally stable at $c/*&o = 0.5. Corresponding flux jump between $/$o ~ and $/$o ~ 1 is 
indicated by vertical line linking the stable solutions. Fig. 4(a). In experiments, the hysteretic loops enclosing these 
lines are usually seen. For tt junctions, Z > Zc, Fig. 4(b), there is only one stable solution with <I>/$o ~ 0.5, and flux 
jumps at ^e/^o = and 1, without hysteresis. For I -C 1, the effect of supercurrents in the loop is negligible, and 
$ « $o for any Z. 

At low temperature, similar results are obtained only in the limiting cases of normal metal and strong ferromagnet 
barriers. For Z = Q and I ^ 1, there is the usual integer flux quantization. Fig. 5(a). For Z = 3 and / ^ 1, there is the 
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half-integer flux quantization, Fig. 5(d), <Ij(<Ijo) being nonhysteretic in contrast to the previous case. Corresponding 
free energies are illustrated in Fig. 6. New features appear for Z close to Z^-, related to the presence of stable and 
metastable minima of the triple-well SQUID energy. Figs. 6(b) and 6(c). For relatively small inductance, I < 10, 
the SQUID has two degenerate configurations: the nearly zero- and one-half- fluxoid state, when ^c/^o < 0.5, and 
the nearly one-half- and one- fluxoid states, when ^c/'I'o > 0.5, as shown in Fig. 7. For example, at T/Tc = 0.1, 
for 0.4 < Z < 1.4 and / < 10, new hysteretic jumps of <& are found, from <I> « to ^/^o « 0.5 in the interval 
< ^c/^o < 0.5, and from ^/^q « 0.5 to ^/^o « 1 in the interval 0.5 < <i?e/$o < 1, Figs. 5(b) and 5(c). This 
anomalous periodicity of $($e) resemble to the usual one. Fig. 5(a), scaled by $o ^o/2. 

IV. TWO-JUNCTION LOOP 

For a loop containing two S/F/S junctions (dc SQUID) the total phase difference is (/> = (/>!— 02, where 0i and 
(/>2 are the phase differences at the first and the second junction, xespectively. In the ground state, the magnetic flux 
through the loop $($c) is calculated from two coupled equationsCJ 

/l($l)+/2($2) =/c (4.1) 

and 

Li/i($i)-i2/2($2)+*-$c=0, (4.2) 

where is the external current, $1^2 = (^o/27r)(/)i.2 and $ = $1 — $2- As expected, $($e) dependence for dc SQUID 
in the absence of the external current, le = 0, is similar to that for rf SQUID. 

At high temperature, T/Tc > 0.9, and for two junctions, Zi < Zc, Z2 < Zc, or two tt junctions, Zi > Zc, Z2 > Z^, 
$($e) dependence is qualitatively the same as in a hysteretic single junction loop. Fig. 4(a). For one junction and 
one TT junction, the two-junction loop behaves as a nonhysteretic single tt junction loop. Fig. 4(b). 

At low temperature, outside the coexistence interval of and tt states, for two junctions, Zi < Z^.- , Z2 < Z^- , 
or for two tt junctions, Zi > Z^+^ Zi > ^c+i the two-junction loop behaves as a hysteretic single junction loop. 
For one junction, and one tt junction, Zi < Z^- , Zi > Z^^ or vice versa, the two-junction loop behaves as a 
nonhysteretic single tt junction loop. For example, at T/Tc — 0.1, $($e) curves for Zi = Zi = and Zi = Zi = 3 
would be practically the same as in Fig. 5(a), and for Zi — and Z2 = 3 as in Fig. 5(d), if Li — Li, the first junction 
normalized inductance li — {2tt/^o)IciLi corresponding to Z of a single-junction loop. When Z, for at least one of 
the junctions, lies within the coexistence interval, Z^- < Z < Z^+ new jumps between the stable solutions appear 
for sufficiently small inductance, similarly to the single-junction loop. For example, at T/Tc = 0.1 and for Li — Li, 
$($0) curves for Zi — and Zi = 0.9 are the same as in Fig. 5(b), and for Zi = and Zi = 1.1, as in Fig. 5(c), see 
Fig. 8. Magnetic flux penetration in the case when metastable states exist in both junctions is illustrated for Zi = 1.1 
and Zi = 0.9 in Fig. 9(a). 

For nonzero external current, up to its maximum value I^ci^e) ^ ^ci + ^c2, 'I'('I'o) curves are similar to those for 
Ic — 0, the flux through the loop becoming closer to the external flux. However, when and tt states coexist, for high 
external current flux jumps appear for any inductance, as illustrated for Ic — Icd^c) in Fig. 9(b). 

In general, /ec(*I'c) curves have characteristic dips at the same positions as the jumps in <l>($c) curves. Three typical 
cases at low temperature Mxe illustrated in Fig. 10 for one period of Icd^c)- For two junction loop, there is the 
usual /cc($e) dependencejl3 with a single dip at ^c/^o — 0.5, Fig. 10(a). For one and one tt junction, there are 
two dips at the ends of the period, <I'c/^'o = and 1, Fig. 10(b). New characteristic behavior is obtained when there 
is the coexistence of and tt states in one of the junctions at least: /ec($e) curve has three maxima and two central 
dips at the flux jumps positions in $($0), Fig. 10(c). 

V. SUMMARY AND DISCUSSION 

The ground state phase difference (/)gs in s-wave S/F/S junctions is always or tt. In the vicinity of the critical 
strength of the ferromagnetic barrier influence Z^, where and tt states are equally stable, the characteristic multinode 
anharmonicity of /((/)) implies the coexistence of stable and metastable tt state for Z^- < Z < Zc, or vice versa for 
Zc < Z < Zc+. For Z < Zc- and for Z > Zc+ only stable states exist, (pgs — 0, and 0gs = tt, respectively. The interval 
of Z, where the metastable states appear, rapidly decreases with increasing temperature, whereas Zc is practically 
unchanged. 
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Stable and metastable states can also coexist in d-wave S/F/S junctions,E3 in contrast to d-wave non-magnetic 
junctions where the similar anharmonicity of the current-phase relation does not lead to the metastable states.EZI 

The peculiar low temperature behavior of clean and transparent S/F/S junctions is due to the modification of 
the Andreev reflection by the exchange interaction in the ferromagnet. The effects of the barrier interfacial non- 
transparency, and of the scattering on impurities, neglected in this paper, suppress the Andreev reflection. As a 
consequence, anharmonicity of-the current-phase relation becomes less pronounced, the value of changes, and the 
metastability region narrows.ca However, when these effects are relatively weak, they do not change the above results 
qualitatively. We also expect that the suppression of the pair potential near theiaterfaces, not included in the present 
non-selfconsistent calculation, would not destroy the metastable states neither .E3 

The coexistence of and tt states in the S/F/S-junction leads to the coexistence of integer and half-integer fluxoid 
configurations in SQUIDs, corresponding to the minima of the triple-well potential energy, and generating two flux 
jumps per one external flux quantum. This is the case when Z is close to Z^ in rf SQUID, and at least in one of the 
junctions in dc SQUID. Another feature is the multipeak dependence of the dc SQUID critical current on the external 
magnetic flux. 

In conclusion, experimental investigation of the predicted new manifestation of the modified Andreev reflection is 
desirable. Beside possible device applications, e.g. quantum computing,c3 the coexistence of integer and half-integer 
flux quantum configuration in S/F/S junction SQUIDs could be also used for the experimental study of the quantum 
superposition of macroscopically distinct states. cJ 
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FIG. 1. Current-phaserelationof S/F/S junction /((/)) for several values of Z, including Zc (dashed curve), for d/^o(0) = 0.25. 
(a) T/Tc = 0.9, Zc = 1.049, and (b) T/Tc = 0.1, Zc = 1.030. 



FIG. 2. Critical current 7c as a function of Z for T/Tc = 0.1 and 0.9. 



FIG. 3. Free energy of S/F/S junction W as a, function of (p for several values of Z, including Zc (dashed curves), for (a) 
T/Tc = 0.9 and (b) T/Tc = 0.1. 

FIG. 4. Magnetic flux $ through the single-junction loop as a function of the external flux $e at high temperatures, 
T/Tc > 0.9, for three values of the normalized inductance I. (a) Z < Zc, and (b) Z > Zc. Only stable (solid curves) and 
metastable (dashed curves) solutions are shown, vertical lines linking the stable solutions. 

FIG. 5. Magnetic flux $ through the single-junction loop as a function of the external flux $e at low temperature, 
T/Tc = 0.1, for three values of the normalized inductance I. (a) Z = 0, (b) Z = 0.9, (c) Z = 1.1, and (d) Z = 3. Only stable 
(solid curves) and metastable (dashed curves) solutions are shown. Vertical lines are drawn to link the stable solutions. 

FIG. 6. Free energy of the single-junction loop VK as a function of $ at low temperature, T/Tc = 0.1, for = 0.1, four 
values of Z, and for <I>e/"l?o ~ (long-dashed curves), 0.5 (solid curves), and 1 (short-dashed curves). 

FIG. 7. Free energy of the single-junction loop as a function of $ at low temperature, T/Tc = 0.1, for = 0.5, and for 
<l?e/<I>o corresponding to the flux jumps at Figs. 5(b) and 5(c): (a) Z = 0.9, ^c/^o = 0.37 (solid curve), $e/*o = 0.63 (dashed 
curve); (b) Z = 1.1, $e/3>o = 0.20 (solid curve), $e/$o = 0.80 (dashed curve). 

FIG. 8. Magnetic flux $ through the two-junction loop as a function of the external flux $e at low temperature, T/Tc = 0.1, 
for Li =1/2, several values of Zi and Z2, and for two values of the first junction normalized inductance h in the absence of 
the external current, le = 0. 

FIG. 9. Magnetic flux $ through the two-junction loop, Zi = 1.1, Z2 = 0.9, Li = L2, as a function of the external flux $e 
at low temperature, T/Tc = 0.1, for two values of the first junction normalized inductance h. (a) le = and (b) le = lec- 

FIG. 10. Maximum external current Joe of the two-junction loop as a function of the external flux $0 at low temperature, 
T/Tc = 0.1, for two values of the first junction normalized inductance h. Three typical cases are shown: (a) Zi = 0, Z2 = 0, 
(b) Zi = 0,Z2 = 3, (c) Zi = 1.1,^2 = 0.9. 
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